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1 Physical setup, assumptions, and notation. 



Consider a flexible pipe, filled with a | fluid under a pressure ~ For the tube not to collapse 



and flow to be possible, it must be | p > p — outside pressure 



we assume 



Po — constant. 



Further, let | x = length coordinate along the tube axis. We now make the following 



Structural Assumptions. 

51. The fluid is incompressible. Let | p = constant be its density. | 

52. The fluid is inviscid (no dissipation by the fluid motion). We also ignore gravity and other 
(possible) body forces on the fluid. The only force acting on the fluid is the pressure. 

53. The walls of the tube/pipe are: 

S3-1. Homogeneous. They have the same properties (thickness, etc.) everywhere. 
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S3-2. Thin and flexible. Thus we will assume that they offer no resistance to bending. Fur- 
thermore, this assumption allows us to ignore the wall thickness, treating it as a surface. 

"Thin" here means this: let L be the scale over which motion occurs — see "dynamic as- 
sumptions" below, and let h be the tube wall thickness. Then L^> h, so that the amount 
of bending in the longitudinal direction is small enough to have negligible effects on the 
force balances. This also has the consequence of keeping the tube cross-section circular — 
see "dynamic assumptions" below, so that there is no bending in the transversal direction. 

S4. The tube perimeter is constant — i.e.: it does not depend on the coordinate x along the 
tube. It then follows (from S3-1) that, under steady state conditions [so that the fluid pressure 
is the same everywhere] the pipe has cylindrical shape. Let 27ra be the repose (no forces) 



perimeter of the tube; i.e. let the repose radius of the tube be a = constant. 



Dynamic Assumptions. 

Dl. Long wave assumption. The | scale L over which the motion occurs] is much bigger 
than the tube diameter — or, for that matter, the tube wall's thickness. This has several 
important consequences: 

Dl-1. Cylindrical geometry. The long wave assumption implies that, at any given point x 
along the tube, the pressure is (essentially) constant — so the cross-section of the tube 
takes a circular shape. Thus, the geometry is locally cylindrical, and can be described by 
the 



radius — a — a(x,t) of the circular cross-section of the fluid region 



inside 



the tube. 



Dl-2. One dimensional approximation. The long wave approximation implies that we can 
neglect any fluid motions in directions transversal to the tube axis. Thus the fluid dy- 
namics can be described by the two scalar functions 
and 



p = p(x, t) — pressure, 



u = u(x,t) — flow velocity along the tube. 



Dl-3. Transversal wall forces only. The long wave approximation implies that the strain 
on the walls is mainly along the perimeter, with relatively little stretching (or bending) 
in the longitudinal direction. Thus, we can assume that the only force by the walls 
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is a "tension" along the perimeter, trying to pull the perimeter back to its 
repose length of 27ra . 



D2. Elastic regime. The amount of wall stretching produced by the motion is small enough 
(and happens slowly enough) that the walls respond elastically. That is to say: the tube walls 
oppose stretching with a force that depends only on the amount of stretching. 

It is important to notice that this assumption involves not just space scales, but time scales 
as well. If the deformations are too large, the walls will not respond elastically — permanent 
deformations will occur, etc. Furthermore, when the deformations occur too fast ( even if they 
are small) dissipation can become important. Of course: what "too large" or "too fast" means 
depends on the physical properties of the tube walls. 

Thus, we can write, for the tension introduced in (Dl-3) a formula of the form: 



where | Aa — a — a | measures the amount of stretching (the tube perimeter changes by 
27rAa), / is a function characterizing the elastic response of the walls, and Tdx is the tension 
force along the perimeter of a "slice" of the tube of length dx. 

D3. Neglect wall inertia. We will assume that the mass of the walls is negligible compared to 
the fluid mass. This will happen if, for example, a^> h, where | h is the wall thickness | — 
this is a stronger assumption than the one in (S3-2), for which only L^> h was required. 

Remark 1.1 Notice that, at least for now, we are not making the assumption of small deviations 
from a steady state 1 . In particular, this would imply that a is nearly constant: a « a s — where 
a s — a is the stretching needed to balance the equilibrium steady state pressure. In this case Hooke's 
law applies, and equation (1.1) can be linearized to 




(1.1) 




(1.2) 



Flow at constant speed and constant pressure. 



Incompressible, inviscid, flow in a narrow elastic pipe with thin walls [R. R. Rosales]. 4 

where 5a — a — a s , | E is the Young's modulus for the wall material, ST = T — T s is the devi- 
ation of the tension from its equilibrium value T s , and \h is the wall thickness"] — which can be 
assumed constant in this approximation. 

Of course, if we do not make the small deviation from a steady state assumption, we cannot assume 
that the wall thickness h is a constant. The variations in h do, of course, affect the forces produced 
by the walls. We can, however, assume that the wall thickness is a function of the stretching, and 
thus we can incorporate the effect of these variations into the force law (1.1) — without being forced 
to track an extra variable h — h(x,t). 

We will make the small deviation from a steady state assumption later. This assumption 
leads to a great simplification of the equations (linearization) . 



2 Derivation of the governing equations. 

We are now ready to derive the governing equations, using conservation of mass and momentum. 
Conservation of Mass. The fluid mass must be conserved. Since S = ira 2 is the cross- 
sectional area of the tube, it follows that the mass density (per unit length dx) is pS, and the 
mass flow is puS. Since there are no mass sources, we must have: 

(pS) t + (puS) x = <=► (S) t + (uS) x = . (2.1) 

By the way: notice that the mass in the tube walls does not "flow". Thus we do not have to 
worry about it when considering the equation for the conservation of mass (or, for that matter, the 
equation for the conservation of momentum). 

Conservation of Momentum. The fluid linear momentum must be conserved. The lin- 
ear momentum density per unit length dx is puS. The linear momentum flux has two 
components: the advective component pu 2 S (momentum carried by the flow), and the mo- 
mentum flux due to the pressure force pS across the tube section. In addition, there is a 
momentum source per unit length, caused by the forces (in the flow direction) on the fluid by 
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the tube walls. It should be clear that this momentum source is given by 2 

M s — 2-Kap—j= x « 2-Kapa x — pS x , (2.2) 
V 1 ~ <>i 

where, from the long wave assumption, we have used that: 1 + a 2 « 1. Thus we have: 

( P uS) t +[pu 2 S + pS) x =pS x <=* (uS) t + (u 2 s) x + ^ Px = 0. (2.3) 

Using (2.1), this last equation can also be written in the form: 

u t + uu x + -p x = . (2.4) 
P 

Equation for the pressure. The pressure in the fluid p = p{x, t) must be balanced by both 
the pressure outside the tube p , and the elastic forces exerted by the tube walls. At the tube wall, 
the pressure force can be decomposed into two components: 

PI. Longitudinal force (per unit area) p -, = = ~pa x . 

P2. Radial force (per unit area) p -, p. 

\J 1 - <>i 

As a consequence of the long wave approximation, the longitudinal force is much smaller than the ra- 



dial one — this is consistent with (Dl-3): the wall elastic forces are transversal only, corresponding 
to stretching mainly in the transversal direction. 

Because the curvature of a circle of radius a is 1/a, the tension T per unit length on the tube walls 
(see equation (1.1)) results in a radial force — per unit area - by the tube walls (see remark 2.1) 
of magnitude 

Forces l -T =±f(^) (2.5) 

per unit area. Because of the assumption (D-3), this force — added to the one caused by the 
external pressure — must balance the radial force (P2) by the fluid in the pipe. Thus, we end up 
with the following formula for the pressure p — p(x, t): 

1 1 /Aa\ , , 

p-p = -T=-f( — ) . 2.6 



2 See (PI): "wall longitudinal force per unit area" x "tube perimeter" = "momentum source" 
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Final Equations. We have now a complete system of equations 



= S t + (uS) 



X 1 



(2.7) 



1 

= Ut + u u x H — p. 

P 



> 



X 5 



where p is given by equation (2.6), and S — ira 2 . 

Remark 2.1 Consider a string under tension, restricted to lie flat in a plane. Let the string be 
described (parametrically) by x — x(s) and y = y(s), where s is the arclength. Let T — T(s) be the 
tension along the string. We now ask the question: What force must be applied to the string, 
to keep it from moving under the effects of the tension? 

In order to answer this question, we first note that the tension produces a force tangent to the string 
at each point, of magnitude T. The tangent vector to the string is given by t — t(s) = (cos#, sin#) 
— where cos(#) = dx/ds, and sin(#) = dy/ds. Thus, the net force on the element of string between 
s and s + ds is given by: 



Now d t/ds — k n, where k — dO /ds is the curvature of the string, and n = (— sin 9, cos 9) is the unit 
normal to the string. Thus the net force on the string produced by the tension can be decomposed 
as follows: 



where the first component is longitudinal, and the second is transversal. This is the force that 
must be balanced to keep the string from moving. In particular, note that if the tension is 
constant, only a force normal to the string is produced — this is the result used above 
to obtain equation (2.5). 

3 Linearized governing equations. 

A simple solution to the governing equations is that corresponding to a steady state, where u — u s , 
a — a Sl S — S Sl and p — p s are all constants. Of course, we must have S s = ir a 2 s and 



Force = T(s + ds) t(s + ds) - T(s) t(s) 



d_ 

ds 



(Tt) ds 



Force = — t + T k n , 

ds 




(3.1) 
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Note that, because the equations are Galilean invariant, we can assume (without loss of generality) 
that u s — 0. 

Consider now solutions that are infinitesimal perturbations of a steady state solution, so that we can 
write: u — u, a — a s + a, S — S s + S, and p s + p — where the variables with tildes are infinitesimal. 
Then 3 we have: 



S — 2 7r a s a and p — — E ha — 



at 



V^Eh ~ _ p 2 ~ 

s~ s cS 



where 




and 



c > is a velocity. Thus the governing equations become 



= S t + S s u x 
c 2 

= ut + —S. 



(3.2) 



where we have dropped the tildes. This yields the 
for both u and S. 



wave equation — u tt — c u xx — S tt — c S x 



THE END. 



3 For the second formula here, see equations (1.2) and (2.5). 



